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Abstract 

We consider a quantum system strongly driven by forces that are periodic in time. The 
theorem concerns the probability P(e) of observing a given energy change e after a number of 
cycles. If the system is thermostated by a (quantum) thermal bath, e is the total amount of 
energy transferred to the bath, while for an isolated system e is the increase in energy of the 
system itself. Then, we show that P(e)/P(— e) = e'^", a parameter-free, model-independent 
relation. 

In the past few years there has been a renewed interest in the study of quantum systems out 
of equilibrium, to a large extent stimulated by the design of new experimental settings and by 
the construction of new devices. If a system is well out of equilibrium, as for example when it is 
strongly driven by periodic forces, then linear response theory (understood as linear perturbations 
around the Gibbs measure) is insufRcient. 

Even in the context of classical mechanics not many generic results are available beyond linear 
response. An interesting new development consists of a number of relations for strongly out of 
equilibrium systems, mainly regarding the distribution of work and entropy production. The first 
of such fluctuation theorems was discovered by Evans et.al who understood that the basic 
ingredient was time-reversal symmetry. 

Two important further steps, made by Gallavotti and Cohen widened the scope and interest of 
the subject. On the one hand, it was realized ||^ that the fluctuation theorems are indeed the far 
from equilibrium generalisations of the well-known equilibrium theorems (fluctuation-dissipation 
and Onsager reciprocity). Most intriguingly, a byproduct of their proof |^ was that, just as the 
validity of the fluctuation-dissipation relation is a strong indication of equilibration, the fact that 
a fluctuation formula holds in a driven stationary system strongly hints that the system can be 
considered 'as ergodic as possible' — with all the implications this entails. 

These results concern deterministic systems. If instead a finite system is in contact with a 
stochastic thermal bath, then the 'ergodicity' questions become trivial, and all the fluctuation 
formulae are extremely simple to prove ^ ^. 
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Particularly relevant for the present work are the simple and remarkable Jarzynski and other 
'work relations ' valid well out of equilibrium [0, ^ |ll| , whose close relation to the fluctuation 
theorems was clarified by Crooks ||l^ . 

Except for Ref.|jll|, the developments described so far are restricted to classical mechanics. 
The purpose of this paper is to prove two versions of a fluctuation theorem for quantum systems 
under strong periodic drive, either isolated or in contact with a thermal bath. Very recently, these 
questions have become relevant in the context of detection of quantum 'shot' noise generated by 
currents flowing through devices, given the possibility it offers to observe different quasiparticle 
charges of the carriers |p^ . 

We shall consider an evolution generated by a time-dependent Hamiltonian H(q,p, t) with real 
matrix elements, so that there is time-reversal symmetry in each infinitesimal time step. We shall 
assume the explicit time-dependence is periodic (with period to) and that the cycles are symmetric 
@: 

H{q,p,t)^H{q,p,t)=H{q,p,-t) (1) 
Splitting the evolution in infinitesimal steps ti, ... t„: 

|(/.(t/)) = Um,)) = UnUn-l ... U2Uim,)) (2) 

where n is large, Ur = e^'^^'' and t = {tf —ti)/n ~> 0, the reality condition implies that: [/^ = U* 
and, together with the symmetry of cycles Eq. (0) this implies that the evolution over a cycle 
satisfies time- reversibility : Consider the case in which the time tf — ti consists of an integer number 
of cycles, and hence Ur ~ C/„_r+i (cfr. (|^)). Then: 

t/t ^ Ulul ... Ul_^Ul = [U1U2 ... U^-iUn]* 

= [[/„[/„_! ... U2U1]* =U* (3) 

a formula only valid for (t/ — ti)/to ~ integer, which wc shall assume throughout this paper. 
We denote the eigenvectors of the initial time Hamiltonian \ipa)'- 

ifll^a) =ea|V'a> (4) 
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and the corresponding partition function: 



Z = Tre"''^^ 

Consider the foUowing protocoh 

• With (canonical) probabihty = e~^^'' /Z we choose a wavefunction 

• We let it evolve through an integer number of cycles: 

10) -L/l^,) =^|V^)(V'^|0) 



(5) 



(6) 



• We measure the final Hi and record the energy difference e between initial and final times. 
Let us calculate the probability distribution of e. For a given value of a, this distribution reads: 

= ^%-(£^-e„)]|(V'7|f^lV^o)l' (7) 

7 

The average distribution P(e) over initial conditions is: 



Pie) = 7 E -^fe - (£7 - e^)] \{num? e-^'- 



(8) 



7.Q: 



Writing the delta function in integral form, we get: 

1 f'ioo 

P{e) = - dAe-^[^-(^-'-^°)lEl(^7|C^IV'o) 

^ J — ioo „, „, 



(9) 



where: 



0(A) = E' 



{lp~f\U\lpa)\ e 



7, a 



E e^^- (^o I C/t I ) (^^ I (7 1 } e- (^+'5) (e„ ) 



= Tr e-'^if/ e 



(10) 



The fluctuation theorem can be proved very simply. Firstly, let us prove the following KMS-like 
[p^ relation. Let D be any real symmetric operator. Then: 



Qd{X) 



Tr 



t/te^^t/e-(^+«^ -Tr er^^+^'^'^Ue^'' 



= Qd{-X* -13)* ^Qni-X-P) 



(11) 



where we have used transposition, cycUc permutation and Eq. (|^). Putting D = iJi, this is a 
form of the fluctuation theorem for Q(A) ~ Qd{^)- 

In order to see what the imphcations of ( ITTI ) are for P{e) we shall need to show that Q{X) is 
analytic on the stripe — /3 < Re(A) < 0. To do this, we use the fact that: 

iTiABf < 2 {TtAA'') (TrSSt) (12) 

Putting A = We^"^ and B = Ue-'^^+l^^"^: 

\Q{X)\' < 2 (Tre-H^^(^)+'^}^^) (xre^R^^^'^^) (13) 

Because the partition function (^) converges for positive temperatures, neither factor diverges if 
-P < Re(A) < 

Inserting (|ll| ) with D — Hi in the analyticity result allows us to shift the integration from 
Re(A) = to Rc(A) = -~/3, and we get: 

P(e) = P(-e)e'^'= (14) 

This is the fluctuation theorem for the probability of energy changes in an isolated system. Note 
that we could have obtained this result directly, without writing the KMS equation, by using the 
(time-reversal) symmetry of U (equation (||)) in (^). 

Taking the expectation value of e~^'' over an integer number of cycles, one obtains a quantum 
version of a Jarzynski work formula [0| : 



= J deP(e)e-^" = J deP{-e) = 1 (15) 

where the overline denotes average over quantum amplitudes and initial conditions. This is a 
rather surprising model-independent result, which we can rewrite as: 

= -ilne^ < -ilne^ = e (16) 

P P 

We obtain the second principle, arising as the familiar inequality 

annealed average < quenched average (17) 
with the initial conditions playing the role of disorder. 



Generalisation and Thermostated systems 



In the most physical setting, we have a system in contact with a bath. The bath can be 
modeled for example with an infinite set of harmonic oscillators coupled to each variable in the 
system. Consider for simplicity the case of one system variable: 

H{x,p^,yi,Pl,...,yM,PM) ^ Hsystern + Hint + Hbath (18) 



with: 



Hbath — ^ 



h muJi — 

2m 2 



H,nt = V C^y^x (20) 



(19) 



and, say: 



2 

H,y,t,„,^ ^ + V{x,t) (21) 



(The interested reader will find an extensive literature on this implementation of heat baths in 
| p5t and references therein.) 

We can define the following bases of eigenvalues: 

HbathlXa) = £l\Xa) 
Osystemli'a') = Oa'\^pa') (22) 

O is any real symmetric operator corresponding to an observable depending exclusively on the 
system variables, with a spectrum bounded from below {e.g. Hsystem{ti)). The wavcfunctions 
\Xa) and li'a) depend only on bath and system variables, respectively. 

We consider an initial condition with no correlations between bath and system constructed as 
follows: 



We choose a bath wavefunction \xct) with canonical probability oc e 



-pel 



We choose a system wavefunction \^/ja') with probability pa' oc e~''°°'. This distribution 
is quite general, given the freedom of choice of O. It could be for example a canonical 
distribution at higher temperature (O = 4-iJi), as resulting from a temperature quench. 



We start with the initial state \xa) ® |'0q') and let it evolve. 



• We measure 



D = H, 



bath 



o 



(23) 



at the beginning and at the end, and record the difference. Note that at this stage the 
operator O which we are measuring is forced by the initial condition we are choosing. 

It is now easy to prove that the fluctuation formula ( p^ holds, with e measuring the difference 
in value of bath energy plus observed value of O. To do this one uses (|ll|) with D as in (|2^). In 
order to guarantee the analyticity of Qd, one must assure that (see argument leading to (^3|)) 

g-2{Re(A)+/3}l5" 



Tr 



= Tr 
X Tr 



g-2{Re(A)+/3}ff6„th 
g-2{Re(A)+/3}o 



< oo 



(24) 



and 



Tr [e^^®''^^-"] = Tr [^^'^'^WHba^th'] q-j. L2Re(A)o" ^ 



(25) 



for —(3 < Re(A) < 0. The traces over the bath variables are bounded, since they correspond to the 
bath partition function. For the initial probabilities of the system, these conditions require that: 



Tr [e-^^'^O] oc Y^P'cT 



(26) 



is finite for < /i < 1, a condition we assume. 

We are now in a position of discussing a general system in contact with a an infinite heat bath 
after many cycles have elapsed (t f — ti oo). If we have that, under these circumstances; 

i) The expected value of O — a property of the system — stays finite as t/ grows. 

ii) The bath receives by virtue of the time-dependent forces an energy proportional to the 
number of cycles ]l6| : e = [tf ~~ ti)eo + ei, with ei finite and and Co independent of the initial 
configuration. 

Then, the derivation above of the fluctuation theorem carries over for long times to any initial 
distribution and to leading order intf — ti, equation (^4|) will be a statement about the probability 
distribution of the energy Co the bath received, i.e. its entropy increase jl^. Measuring the 



fluctuation formula in a concrete situation becomes then a test for a property of the 'stationary' 
(i.e. periodical) asymptotic quantum state. 
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